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Abstract. Phyllosilicate is a sheet of silicate tetrahedra bound by 
basal oxygens. A phyllosilicate excitable automaton is a regular net- 
work of finite state machines, which mimics structure of a silicate sheet. 
A node of the silicate sheet is an automaton, which takes resting, excited 
and refractory states, and updates its state in discrete time depending 
on a sum of excited states of its three (silicon automata) or six (oxygen 
automata) closest neighbours. Oscillator is a localised compact con- 
figuration of non-quiescent states which undergoes finite growth and 
modification but returns to its original state in a finite number of steps. 
We show that phyllosilicate excitable automata exhibit spiral and target 
waves, and oscillating localisation dynamics. Basic types of oscillators 
are classified and characterised. 
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1. Introduction 

Phyllosilicates are parallel sheets of silicate tetrahedra, they are widely 
present in nature and typically found in clay-related minerals on the Earth 
surface (TDJIHIS]. Phyllosilicates are used in development of nano-materials, 
nano- wires and patterned surfaces for nano-biological interfaces [UJ EJ [19] . 
They are also employed, in a form of cation-exchanged sheet silicates, as cat- 
alysts in chemical reactions [2], e.g. nickel phyllosilicate catalysts |TT | [T3 |, fT8j . 
In the paper we continue lines of enquiry into space-time dynamics of cellu- 
lar automata on non-orthogonal and aperiodic lattices, including triangular 
tessellations and Penrose tilings [31 [151 [8]. We define an automaton network, 
where connections between finite automata are inspired by simplified struc- 
ture of silicate sheets (as lattices of connected tetrahedra), and investigate 
the dynamics of excitation on the automaton networks for various excitation 
rules. 

The automata models studied in the paper are abstractions of silicon 
sheets. The automata models do not aim to compete with mainstream com- 
putational chemistry models [17]. What are advantages of the automata 
approach? Automata models is a fast prototyping tool for an express eval- 
uation of a space-time dynamics of a spatially-extended active nonlinear 
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Figure 1. Phillosilicate automata. Automata representing 
silicon are small discs, oxygen automata are large discs. Un- 
shared oxygen atoms of the tetrahedra are not shown. 



medium for different excitation rules, and prototyping of unconventional 
computing devices based on the nonlinear medium. Even abstract three- 
state automata are not totally pointless abstractions in imitating excitations 
in active media. Consider, e.g., Belousov-Zhabotinsky (BZ) medium. All 
types of travelling waves — spiral and target waves in excitable mode and 
wave-fragments (dissipative solitons) in sub-excitable mode are sufficiently 
imitated in Greenberg-Hasting automata [9]. The automaton model of BZ 
medium is phenomenologically equivalent to more accurate, close to physical 
and chemical reality, Oregonator model [7] . By analogy, studies of phyllosil- 
icate automata may help us to get an insight into dynamics of localizations, 
defects and excitations, in the lattices of silicon tetrahedra. By developing 
automata models of excitation in tetrahedra sheets we contribute towards 
future developments of silicate sheet based computing circuits where quanta 
of information are transferred and processed by stationary and travelling 
localised excitations pQ. 



2. Phillosilicate automata 

In silicates, coordinated tetrahedra of oxygen anions are formed around 
each silicon cation. In late 1920s Pauling proposed that phyllosilicates are 
sheets of coordinated (SiC^) 4- tetrahedra units, where each tetrahedron 
shares its three corner (basal) oxygens of neighbouring tetrahedra [161 E21 E] 
(Fig.[l]). Vacant oxygens of all tetrahedra point outside the lattice, away of 
tetrahedra; these apical oxygens are not taken into account in our present 
simulation. 

A phyllosilicate automaton A is a two-dimensional regular network Fig. [I] 
of finite state machines. There are two types of automata in the network: 
silicon automata s (centre vertex of each tetrahedron) and oxygen automata 
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(e) a\s) > 1, <t\o) > (f) a\s) > 1, a\o) = 1 



Figure 2. Configurations of excitation developed after ex- 
citing a single oxygen node in the resting lattice. Predicates 
C x (a t (x)) are show in sub- figures' legends. Excited nodes are 
solid discs, refractory nodes are grey discs, resting nodes are 
empty discs. Configurations are recorded at 25th step of the 
lattice iteration. 
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o (corner vertices of each tetrahedron). A silicon automaton has three neigh- 
bours. An oxygen automaton has six neighbours (Fig. [I]). The automata 
take three states • (resting), + (excited) and — (refractory), and update 
their states s f and o l simultaneously and in discrete time t depending on 
states of their neighbours. Let a f (x) be a number of excited neighbours of 
automaton x, < cr*(s) < 3 and < cr f (o) < 6. Each automaton- node 
x — 5, o of A updates its state x l by the following rule: 



(1) 



, if x l = • and C x (a t (x)) 

, if x l = + 

otherwise 



where predicate C x (a\x)) e {(^(x) > 0), (a\x) = 1), (a\x) > 1), (a\x) = 
2)}- 

Exemplar configurations of excitation developed after single-node excita- 
tions are shown in Fig. [2j A typical 'circular' wave of excitation develops 
when C(sY = (a\s) > 0) and C(o) 1 = (a\o) > 0). (Fig. §i). Pat- 
terns of target waves are observed in situations C(sY = (^(s) = 1) and 
C{of = (a\o) > 0) (Fig. |^), and C(sY = > 1) and C{of = 

(cr^(o) > 0) (Fig. [2^). Wave-fragments are formed when C(o)* = (cr^o) = 1) 
and C(sY e {(a\s) > 0), (a\s) = 1), (a\s) > 1)} (Fig. §bdf). 

In scenario illustrated in Fig. [2)3 we observe fragments of excitation waves, 
open ended segments of wave-fronts. These wave-fragments expand. Their 
ends bend backwards, split away, new wave- fragments propagate towards 
originally excited site. They collide with outcome wave- fragments, new 
wave-fragments are produced in the result of collision. Being strictly de- 
terministic an excitation dynamics inside expanding perturbation domain 
might look chaotic (Fig. [2)3). In (Fig. [2jl) wave- fragments undergo a bi- 
nary division at every iteration. This is because only resting nodes with 
exactly one excited neighbour are excited, and therefore only end of any 
wave-fragment transfer excitation to the resting medium; interior of the 
wave- fragment extinguishes. 

Development of A after random (each node gets excited with probability 
0.3) excitation for all possible combinations of C f (x), x = o, s is classified in 
Tab. [TJ For excitation rules cr f (o) > and 0^(0) = 1, a random excitation 
leads to formation of generators of waves and wave- fragments. For ^(s) > 1 
and o~*(o) > 1, and o~*(s) > 2 or o~*(o) > 3, all excitation, evoked by a random 
stimulation, extinguishes after few iterations and the automaton returns to 
its resting state. 



3. Oscillators 

For excitation functions o~*(s) > and 0^(0) > 1 (/01), ^(s) > and 
cr\o) = 2 (/ 02 ), a\s) = 1 and a\o) > 1 (/ n ), and a\s) = 1 and cr\o) = 2 
(/12); a randomly excited automaton A evolves towards a configuration with 
few localised excitations. 



ON OSCILLATORS IN PHYLLOSILICATE EXCITABLE AUTOMATA 



5 



Table 1. Outcomes of random excitation of A. Excitation 
does not persist in scenarios marked by 
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Figure 3. Average parameters of oscillators, (a) Distribu- 
tion of maximum linear sizes, (b) Distribution of numbers 
of excited states in oscillators. The distributions were calcu- 
lated on 1677 oscillating patterns (not necessarily different 
types of oscillators) discovered in 40000 trials. 
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(m) (n) 



Figure 4. Exemplar configurations of oscillators with a to- 
tal, i.e. over all configurations of an oscillator, number of 
excited states 12 (a-f), 14 (g-j) and 16 (k-n). All oscillators 
shown have period three, their configuration change from left 
to right. Excited state is black with red contour, refractory 
state is grey/blue. 
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In terms of cellular automata, these excitations are oscillators — compact 
configurations of non-resting states undergoing modifications yet returning 
to their original states in a finite number of steps (Tab. [I]). 

To discover how often oscillators appear in randomly excited automaton 
we stimulated A with a rectangular domain of 10 by 10 nodes. At the 
beginning of evolution each node in the domain is assigned one of three 
states equiprobably. In each trial, the automaton was iterated until all 
excitation disappeared or a localised excitation emerged. For each function 
foil /o2? fn an d fi2 10K trials were conducted. We found that oscillators 
appear with probability 0.057 in automata governed by function /qi, 0.086 
in automata governed by /n, 0.0144 for /02, and 0.0108 for fi 2 . 

Distributions of oscillators on linear sizes and weights (numbers of excited 
states) are shown in Fig. |3j The distributions are obtained by averaging the 
maximal linear size and the number of excited states over all the configu- 
rations of a period of the oscillator. A typical oscillator spans five to eight 
nodes along one dimension and consists of four to seven excited states. Ex- 
amples of minimal size and weight oscillators are shown in Fig. [4j All oscil- 
lators discovered exhibit the same behaviour for all oscillator functions. All 
oscillators studied have period three. Based on a total number e of excited 
states, summed over three configurations of an oscillator, and a sequence of 
excited states w = (w t w t+1 w t+2 ) , where w t is a number of excited nodes 
at step t, in oscillator cycle, we can select the following classes of minimal 
oscillators: 

• Class C e= i2. These are lightest oscillators, see examples in Fig. [I| 
There are three sub-classes determined by w: w\ 2 = (444) (Fig.[4p,d), 
w \ 2 = (345) (Fig. gbc) and w\ 2 = (354) (Fig. gjjf ) 

• Class C e= i4. These are middle weight oscillators. There are three 
sub-classes: w\ A = (446) (Fig. [4^), and w\ A = (455) (Fig. [ZJii). 

• Class C e= iQ. These are heavy weight oscillators. There are two 
sub-classes: w{ 6 = (466) (Fig. [Z|n) and w\ Q = (556) (Fig. [IJan). 

Ratio of excited silicon automata to excited oxygen automata in oscilla- 
tors configurations are 0:3, 2:2 and 4:1 in class C e= i2] 1:3, 3:3 and 5:1 of 
class C e =u; and, 0:4, 2:3 and 4;4 in class C e =iQ. In class C e =i2 excitation 
is equally distributed between silicon and oxygen automata, due to equal 
chances of oscillator configurations with ratios 0:3 and 4:1 to appear in a 
randomly perturbed automaton. Silicon automata are excited more often in 
oscillators of class C e= i4, assuming equiprobability of the oscillator genera- 
tion. Oxygen automata are more likely to be excited in oscillators of class 

C e =16- 

We have analysed in details only oscillating localizations with minimal 
size and number of non-resting node-states. In principle, the only limit on 
oscillator size is a size of automaton network. Examples of large oscillators 
are shown in Fig. [5j Oscillator hi change its maximum linear size between 
eight and nine, and consists of five, six and nine nodes in each of its three 
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(a) hi 



(b) h 2 




(c) h 3 



Figure 5. Examples of heavy oscillators. Excited state is 
black with red contour, refractory state is grey/blue. Each 
oscillator has three configurations, shown in their succession 
from left to right. 



configurations (Fig.[5^). Silicon and oxygen nodes are excited in total equal 
ratio. Oscillator I12 spans over 13 nodes in its largest configuration. Each 
configuration of 1i2 consists of five excited silicon nodes and five or six excited 
oxygen nodes (Fig. [5)3). Oscillator /13 is the largest and heaviest amongst 
examples in Fig. [5j its maximum size is 17 nodes and it has 15 and 17 
excited sites in its configurations. Oxygen automata are slightly more often 
excited than silicon automata in patterns of /13; exact ratios of s to o are 
6:9, 8:9 and 9:8. 

Oscillators survive change of excitation rule from an oscillator supporting 
rule to a wave supporting rule and back to an oscillator supporting rules. 
When a node state transition rule is dynamically changed to a wave sup- 
porting rule an oscillators becomes a source of waves generation (Fig. |6]l). 
Resting nodes adjacent to excited sites of an oscillator becomes excited and 
excite their neighbours at the next time step. Thus waves of excitation start 
to emerge around the original oscillators. The wave propagate outwards 
(Fig. [6^-n). When the node state transition rule is switched back to the os- 
cillator supporting rule (Fig. ^jp) the excitation waves extinguish (Fig. |6]3-q) 
yet original oscillator remains active (Fig. (6t-t). 
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(m) t = 22 (n) t = 23 (o) t = 24 (p) t = 25 

Figure 6. Dynamics of automaton A, where excitation rules 
are changed twice during the automaton's development. Au- 
tomaton evolves by rule /oi for t = 1, . . . , 4. At step it — 4 
excitation condition is modified as cr*(o) > (to make the 
wave supporting rule). At step t — 20 the rule is changed 
back to oscillator rule /oi- Excited states are black with red 
contour, refractory states are grey/blue. 
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4. Discussion 

We modelled a network of excitable finite state machines based on a 
simplified topological structure of silicate tetrahedra sheets. We studied 
a domain of excitation rules that support oscillating localised excitations, 
provided basic classification and analysed structure of minimal oscillators. 
The following topics remain open. 

All oscillators discovered so far have period three steps. Are there oscil- 
lators with period more than three? 

We did not find gliders, travelling localizations, which would propagate 
for a long distance without splitting into new localizations. Is it because the 
rule does not support gliders at all or we were unlucky in not discovering 
a glider in 100K random perturbations of phyllosilicate automaton? 
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